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Abstract 

The standard semiclassical calculation of transmission correlation functions 
for chaotic systems is severely influenced by unitarity problems. We show 
that unitarity alone imposes a set of relationships between cross sections cor- 
relation functions which go beyond the diagonal approximation. When these 
relationships are properly used to supplement the semiclassical scheme we 
obtain transmission correlation functions in full agreement with the exact 
statistical theory and the experiment. Our approach also provides a novel 
prediction for the transmission correlations in the case where time reversal 
symmetry is present. 
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I. INTRODUCTION 



Since the pioneering work of Bliimel and Smilansky |T| the semiclassical S-matrix 
has been used by many authors to study fundamental questions related to quantum chaotic 
scattering. In recent years the interest in this subject has grown due to the experimental 
investigation of electronic transport through small devices, such as open quantum dots [§. 
At sufficiently low temperatures, these devices preserve quantum coherence and are called 
mesoscopic. The conductance in transport processes that preserve quantum coherence is 
directly related to the S'-matrix by the Landauer-Biittiker formula [^. Since the classical 
underlying electronic dynamics in quantum dots is believed to be chaotic, these are excellent 
systems to observe the quantum manifestations of classical chaotic scattering 

One of the central issues in mesoscopic physics is to single out statistical universal prop- 
erties of quantum systems. In this paper we show that this goal cannot be theoretically 
achieved using the standard semiclassical S-matrix approach, since the latter is not able to 
provide trustworthy results for universal cross-section or conductance correlation functions. 
We also show that this situation can be fixed by imposing a set of semiclassical sum rules 
that guarantee unitarity. 

For a device connected to reservoirs by two leads, the Landauer-Biittiker formula relates 
its conductance G to the transmission coefficient T by the expression G = {2e'^/h)T. In elec- 
tronic transport, T is usually called the dimensionless conductance. When the entrance and 
exit leads support A^^i and N2 propagating modes or channels, respectively, the transmission 
T reads 

Ni N 

T{x) = Y: E \Sa,{x)\' = Y.^a,{x), (1) 
a=l b=Ni+l a,b 

where we have introduced aab, the transition probability (the cross- sect ion, apart from a 
kinematical factor) from channel b in one lead to channel a in the other lead. The parameter x 
represents either the energy E, or an external parameter X such as a magnetic field, or both. 
Without loss of generality we restrict our discussion to the case where Ni = N2 = N/2 = M. 
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The original problem of conductance fluctuations in a mesoscopic device is now cast into 
a more generic one of transmission fluctuations of a quantum process where a particle is 
chaotically scattered. Thus, this discussion is also of interest for transmission experiments 
with irregular microwave cavities, "chaotic" atoms and nuclei. Our goal is to use the semi- 
classical theory to describe the statistical properties of the transmission as the parameter 
X is varied. This information is contained in the average value of T (a two-point statis- 
tical measure of the S'-matrix elements) and in its autocorrelation function (a four-point 
function), defined as 

{TX)^{nx)T{x'))-{T{x)){T{x')). (2) 

The average is taken over x and x' keeping the difference |a: — x'| fixed. The transmission 
autocorrelation function is directly related to the covariances of the transition probabilities 

M N 

iT,T')= E E i^at,a'J. (3) 

a,c=l b,d=M+l 

The variance of T, var(T), is the statistical measure of a fundamental phenomenon in 
mesoscopic physics: For systems where quantum coherence is preserved the dimensionless 
conductance displays fluctuations of order unity irrespective of sample size, provided the 
dynamics is chaotic (or diffusive) and there are no tunneling barriers hindering the trans- 
mission. This phenomenon is known as "universal conductance fluctuations" (UCF) [Q. 
Thus, a successful approximation scheme to explain UCF has to be accurate to the level of 
unity for the variance of T. In the specific case of quantum dots, i.e., ballistic electronic 
cavities, the random matrix theory |]^,|^ and the supersymmetric method are, so far, the 
successful approaches to calculate var(T). 

The purpose of this work is twofold. First we analyze a very simple statistical measure - 
the average cross section - to show that the standard semiclassical S'-matrix theory does not 
achieve the required precision to be useful as a theory for UCF. In doing this, we indicate its 
main sources of inaccuracy and discuss the main problems involved in improving the theory. 
We then show how to fix the inaccuracies by making explicit use of the unitarity of the 5" 
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matrix. This procedure is similar in spirit to those used in semiclassical studies of spectra of 
closed system [|10|. It can be viewed as a proposal for a set of semiclassical scattering sum 
rules to enforce well known exact symmetries of the S'-matrix. 



II. THE SEMICLASSICAL APPROACH 

We start with Miller's semiclassical S'-matrix formula (including now transmission 
and reflections) 

S,,{E,X) = Y: ^Jp,{E,X)e^'^-^^'''^/^ , (4) 

/i(a,f)) 

where fi{a, b) labels the classical trajectories that start at channel b and end at channel a, 
(j)fj_ are their reduced actions (with a Maslov phase included), and stands for the classical 
transition probability of going from a to 6 through the orbit fi |Tl[] (here and throughout the 
paper the tilde indicates that the semiclassical approximation is employed). It is implicit 
in the derivation of Eq. (|^) that the number of open channels must be much larger than 
one and that there are no tunneling barriers between the leads and the cavity. When the 
scattering is chaotic (and the short time dynamics does not significantly contribute to S) the 
domains of applicability of the semiclassical theory and random matrix theory coincide and 
both approaches should be comparable. In this regime the semiclassical approach provides 
the dynamical explanation for the universality of the scattering fluctuations. 

In general the semiclassical S-matrix is not exactly unitary at any given energy E. 
Indeed, it is only upon energy averaging and for the case of broken time-reversal symmetry 
(BTRS), that unitarity is automatically fulfilled. We now show this known result so as to 
present the basic elements and approximations employed in this paper. The energy averaged 
semiclassical cross section reads 

{aat) = {\Sat\') = E (v^e^^*"-'^''^/') • (5) 

Here (■ ■ ■) indicates an energy average within an energy window where the classical dynamics 
presents little changes, nonetheless comprising many resonances. To compute the energy 
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average one neglects the energy dependence of the probabihties and uses the diagonal 
approximation. The latter says that, on average, only orbits having the same action are 
correlated. If there are no symmetries present, this means that (exp[i(0^ — = 6^u- 

Then 



The proof is completed by using the classical normalization condition ||TT 

N 

E E P. = 1 . (7) 

a=l ,i{a,b) 

which insures that J2a{\^ab\^) = 1- 

We shall assume that for any given entrance channel b, all exit channels a are equivalent, 

i.e., 

E P. = ^ , (8) 

/i(a,6) 

which yields 

i^ab) = ^ . (9) 

The assumption of equivalent channels is justified (in the BTRS case) if the particle typically 
stays inside the interaction region time enough to be randomized, meaning that it becomes 
equiprobable to be ejected through any outgoing channel. The analysis of the other limiting 
case where time reversal symmetry is preserved will be postponed to Section |T|. 

In order to clearly explain why unitarity problems affect the semiclassical theory of 
transmission fluctuations we introduce the object 

_ N N 

la = E = E E V^e^^^-'^")/^ . (10) 

6=1 6=1 

f(a,t») 

If the semiclassical S'-matrix had been exact 1^ = 1; instead one has = 1. The lack 
of precision of the standard semiclassical scattering theory at the four-point level becomes 
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evident by analyzing the variance (la, Ife). The semiclassical approximation gives (Iq, 1;,) 7^ 0, 
leading to a "unit fluctuation" problem. To see this, using Eq. (H), we write 

(la. If.) = ( H H H VPtJ^^W^ exp [U(pf, -(j)u + 0m' - 0^')] ) - 1 (11) 

i'{a.c) v'{b,d) 

SO that the diagonal approximation yields 

N 

(laAb) = J2 J2 + '^M'^"^!^/^') - 1 • (12) 

C,d=l fJ,(a,c) ^i'(b,d) 

The first Kronecker 6 product decouples the sums over orbits starting at channel c and 
ending at a from those entering the scattering region through channel d and exiting through 
b. The resulting double sum adds up to unity. The second product contains crossed terms 



which vanish unless a = b and c = d. Equation (JT^) becomes 



^ 6 

(la, h) =5abJ2 J2 P^^P^ = ■ (13) 

v(a,c) 

This inaccuracy is neither unexpected, nor large. However, it has important consequences 
for the calculation of transmission fluctuations. This becomes evident by inspecting 

/ N ^ N ^ \ / N N \ 

\a=l 6=1 / \a,c=l b,d=l J 

which has the same double sum structure of the transmission variance and shows an inac- 
curacy exactly of the order of the effect that we aim to describe. 

Let us be more explicit and go back to the analysis of the transmission autocorrelation 
function. Recalling Eq. and assuming channels to be statistically equivalent we write 

(T, T') = M'iaat, cr'J + 2M\M - l)(a,,, <J + M\M - ifia,,, a' J . (15) 

Here we use the convention that different indices a and b in the covariances imply that a^b. 
This means that for the above equation 6 7^ c in the second term of its RHS, and a ^ c 
and b ^ d for the third one. We shall demonstrate below that, owing to unitarity, all three 
terms in the RHS of Eq. ([T5|) are of the same order of magnitude. However, within the 



diagonal approximation, both the semiclassical covariances icfab,d''^c) i^ab, o'cd) are zero. 
Let us admit that the semiclassical approach gives the correct result for the transmission 
probability covariances to order 1/N'^. Then, for a successful description of the transmission 
fluctuations, the theory has to be improved to access the first non vanishing order in the 
non diagonal terms {crab,o'ac) {crab,o'cd)y which are 0{N~^) and 0{N~^), respectively. 
Though desirable this is not really necessary. The alternative scheme we propose is to bypass 
the explicit semiclassical calculation of the nondiagonal covariances and use the unitarity of 
the S matrix to relate the latter covariances to the diagonal one. Having expressed (T, T') 
in terms of {aab, o-'^b) alone, we use the semiclassical approximation only at the very end. 



III. ENFORCING UNITARITY 

The relations among diagonal and nondiagonal covariances can be easily obtained from 

N 

{J2 ^'^b, (y'cd) = (i> (^cd) = , (16) 

6=1 

which follows from unitarity. To proceed further we have to separately analyze the cases 
where either time reversal symmetry is absent (BTRS) or present (TRS). This distinction 
is necessary because the "elastic" processes (a = b) and the "inelastic" ones (a 7^ b) display 
different statistical properties when time reversal symmetry is preserved. Indeed it is well 
known that due to quantum interference, time reversal symmetry enhances the average 
reflection probability (aaa) by a factor of two |T^. This can be understood semiclassically 



by noting that in the TRS case there are pairs of orbits having the same action (time reversal 
partners) which contribute to aaa, thus interfering constructively to produce the factor two. 
Due to this effect the classical equivalence of channels breaks down at the quantum level. 
In this case the equivalence is only restored when time reversal symmetry is broken, leading 
to the considerations presented in Section 0. 

To illustrate the consequences of this phenomenon and the spirit of our scheme, let us 
analyze {T{x)) in the crossover regime from TRS to BTRS. Recall that 
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M N 

{T{x)) = Y: E {^a,{x)) = M'{a^,{x)). (17) 

a=l b=M+l 

Here x parameterizes a Hamiltonian change breaking time reversal symmetry as x grows from 
zero (TRS) to some critical value x* (BTRS). Unitarity relates diagonal and off-diagonal 
averages: 

l = {aaaix)) + iN-l){aabix)) . (18) 

This equation allows us to write (T) in terms of {caa), the average that is semiclassically 
sensitive to time reversal effects, to obtain 

{Tix)) = ^^^^(1 - i^aaix))) . (19) 

For X = the elastic enhancement is maximal and hence (T(x)) takes its smallest value. 
In mesoscopic physics, to distinguish from strong localization which is a phenomenon very 
different in origin, the reduction of transmission due to TRS is called the weak localization 
peak. Up to this point Eq. ([T9| ) is an exact expression. The semiclassical result is obtained 
by calculating {craa{x)) from Miller's formula. For x representing a magnetic field, the semi- 
classical approach gives a Lorentzian shape for the weak localization peak [|13[ in agreement 
with RMT ||14[. The amplitude of the weak localization correction can be readily obtained 
recalling that (craa(O)) ~ 2/N and {aaa{x > x*)) ^ so that 

(f (0)) - (f (x*)) = -\ , (20) 

again in agreement with random matrix theory [0. The discussion above is not entirely 
original and was inspired by the pioneer semiclassical study of the weak localization peak in 
ballistic cavities developed by Baranger and collaborators [|13|. Based on the same strategy 
presented above, we are now ready to understand the UCF problem. 

A. Transmission fluctuations in systems with broken time reversal symmetry 

In this case all S matrix elements are statistically equivalent. In order to express (T, T') 
in terms of {aat, cr'^b) it suffices to consider the following two independent unitarity equations 
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TV N 

<y'ac) = , ^{<yab, <y'cd) = 

6=1 b=l 

(c 7^ a in the second equation). The above relations can be reduced to: 



(21) 



(^afe, o'J + {N- 1)((T,,, a' J = , {aab, a' J + (iV - l)(a,,, a' J = 



(22) 



(the convention about indices being the same as in Eq. (|15D). Notice that these relations 
are not satisfied in the diagonal approximation. At the semiclassical level Eqs. (^Tj) and 
(^) can be regarded as sum rules that go beyond the diagonal approximation. Insertion of 
Eqs. ([2^ ) into Eq. (pIS]) readily renders 



(T,T') 



(23) 



(2M- 1)2 

As in the two-point analysis this equation is exact. 

Now we are ready to employ the semiclassical approximation to compute (o"a{,, a^^). Let 
us first consider the case where x stands for the energy, i.e., 



[(Tab, Cr'^b) = Cab{e) = { aab{E + |) (TabiE - |))^ " {^ab)l ■ 



(24) 



The semiclassical autocorrelation function Cabi^) can be calculated for classically small val- 
ues of e, i.e., for energy differences such that the classical perturbation theory holds. In this 
case one keeps the stability coefficients constant and expands the actions to first order in 
e, i.e., (f)fj,{E ± e/2) ^ 4>fi{E) ± here is the time the particle takes to travel from 

channel b to channel a along the orbit fi. After the diagonal approximation we obtain 



1^(0,6) v'(a,b) 



X exp 



le 



2^ yTf^-T^ + T^' - V) 



(25) 



Using the same arguments employed after Eq. (|^) we arrive at 



fi{a,b) 



(26) 



According to the analogue of the Hannay-Ozorio de Almeida sum rule for open systems [jTS 



E P. = ^e-^*5t, (27) 

t<T^<t+St 

where J2t<T^,<t+5tPiJ. is the sum of all classical transition probabilities from channel a to 6 
through trajectories within a small time interval [t,t + 6t], where 6t is classically small. The 
exponential is determined by the inverse escape time 7, later to be associated to an energy 
width r = h'j. Replacing the sum over orbits by an integral over the time, we finally obtain 

11 

Ca>.ie) = YTWy ■ ^^^^ 

Before commenting on this result, we shall generalize it by accounting for an external 
parametric change X in the Hamiltonian. For instance, in (mesoscopic physics) experiments 
X is frequently an external magnetic field. The strategy to compute the generalized cor- 
relation function is, as above, based on classical perturbation theory. By expanding the 
reduced action to first order in X, we have to deal with = d(j)^/dX. A full account 
of the technical details involved in calculating the parametric correlation can be found in 



Ref. [|I6|. The basic step though is to compute the time average 

(29) 



e*0-'5^/";^^=exp 



2r 

Since {Q'^{t))st grows diffusively with time, i.e., {Q'^{t))st = o^t, the autocorrelation function 
becomes 

-11 

Cabis, 5X) = j;p^^^^^x/X,)^]^ + {e/Ty ' 

with 

X^ = 2hV/a . (31) 
Then the semiclassical result that includes unitarity restrictions is 
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Remarkably, this result agrees exactly with the dimensionless conductance autocorrelation 
function for open ballistic dots in the limit ^ 1 obtained by Efetov [Q using the super- 
symmetric technique. The agreement extends also to the structure of the parameter Xc (|3lD 
if one relates a to the level velocity of closed systems as defined in Ref. [0. In this respect 
the semiclassical approach is complementary to random matrix theories, since it provides a 
dynamical interpretation for the nonuniversal quantities Xc and F. 



B. Systems with time reversal symmetry 

Now the S matrix is symmetric and the statistical properties of diagonal and off-diagonal 
elements are different. Hence we need to write down two additional unitarity relations, as 
compared with the BTRS case, in order to single out the elastic case separately: 

N N N N 

II(t^a6, Cr'^a) = , Y.(^cb, CT^J = , ^(cTafe, CX^J = , JK^^b' ^cd) = (33) 
b=l 6=1 6=1 6=1 

(c, d ^ a). In terms of the basic covariances, the above system of equations is rewritten as 

(f^aa, (y'ab) + (^aa, (^'J + {N - 2){aaa, Cr[^) = 
{(^ab, (Tab) + {(^aa, (^ab) + {N - 2){aab, (T^J = 

{a^a, ai) + 2{aab, a'J + (iV - 3)(a,fe, a'J = , (34) 



where the channel index convention is that following Eq. (|15|) . The only difference with the 
BTRS case is the factor 2 in the last equation, which is a consequence of the symmetry of S. 
As in the preceding subsection, we would like to express all covariances in terms of {(Jaa, cr'aa) 
and {aab, (y'ab) which are the only nonzero ones in the diagonal approximation. Regrettably, 
we have more unknowns than equations, and it is not possible to obtain an exact equation 
like Eq. (^). However, as we are only interested in a relation which is correct to leading 
order in it suffices to consider the simplified system 

{aab, o'J + N{aab, a'J = 0{N-') , 2(a,,, a' J + Ar(a,,, a' J = 0{N-') , (35) 
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which is obtained from the last two equations in (|3^ by keeping only the leading terms. 
These relations are the TRS analogues of ( ^2]) and lead to 

(T,r)^^K,,<j. (36) 

Semiclassically, time reversal effects only manifest themselves in the diagonal covariances 
{aaa,o''aa)- The correlator {(Tab,cr'^b) is the same as that for the BTRS case. Thus, at the 
semiclassical level, the effect of time reversal symmetry is to enhance the transmission fluc- 
tuations by a factor of two (cf. Eq. (^)), without changing the shape of the correlation 
function, that is 

f ') = IjrusxJxW^Wy ^ "^""^ ■ 

We are not aware of any study of the transmission autocorrelation function for the TRS 
case using the supersymmetric method. For technical reasons it is much simpler to obtain 
the variance var(T) using random matrix theory. This result is known and, for ^ 1, 
agrees with our semiclassical calculation 

var(f ) = - . (38) 



IV. CONCLUDING REMARKS 

Spectral studies in closed system have shown that the diagonal approximation should 
start to fail when the orbits involved have periods of the order of the Heisenberg time th- 
In scattering systems, the contribution of orbits with periods larger than the mean escape 
time Te is negligible. Given that in the semiclassical regime Tg <C th it is generally accepted 
that the diagonal approximation should be unproblematic for scattering systems. 

However, we have shown that the standard semiclassical approach fails to describe the 
transmission fluctuations because the diagonal approximation does not preserve the unitarity 
of the S matrix to the required precision. One way to circumvent this problem, perhaps 
the most satisfactory from a theoretical point of view, is to improve the semiclassical theory 
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to include correlations between different orbits. Alternatively, we have shown that the 
unitarity of the S matrix can be used to express the transmission autocorrelation function 
in terms of transmission probabilities. Such expression contains some information about 
unitarity allowing the standard semiclassical approximation to be invoked resulting in a 
theory consistent with UCF. 

Other difficulties are also encountered when calculating the transmission correlations: 
The semiclassical transmission correlator is not translationally invariant. For instance, this 
is manifest in the fact that 

{f{E)f{E + e)) ^ {f{E + e/2)f{E - e/2)) , (39) 

which can easily be checked by inspection. In our calculations we preferred to use {T{E + 
e/2)T{E — e/2)) because it is exphcitly real. This choice is consistent with the spirit of 
this work, i.e., all information about exact quantum symmetries must be used in trying to 
compensate for the shortcomings of the semiclassical S matrix. 
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